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Abstract: Levodopa (LD) is used to increase dopamine level for treating Parkinson's 
disease. The major metabolism of LD to produce dopamine is decarboxylation. In order to 
understand the metabolism of LD; the electronic structure of levodopa was investigated at 
the Density Functional DFT/B3LYP level of theory using the 6-311+G** basis set, in the 
gas phase and in solution. LD is not planar, with the amino acid side chain acting as a free 
rotator around several single bonds. The potential energy surface is broad and flat. Full 
geometry optimization enabled locating and identifying the global minimum on this 
Potential energy surface (PES). All possible protonation/deprotonation forms of LD were 
examined and analyzed. Protonation/deprotonation is local in nature, i.e., is not transmitted 
through the molecular framework. The isogyric protonation/deprotonation reactions seem 
to involve two subsequent steps: First, deprotonation, then rearrangement to form 
H-bonded structures, which is the origin of the extra stability of the deprotonated forms. 
Natural bond orbital (NBO) analysis of LD and its deprotonated forms reveals detailed 
information of bonding characteristics and interactions across the molecular framework. 
The effect of deprotonation on the donor-acceptor interaction across the molecular 
framework and within the two subsystems has also been examined. Attempts to mimic the 
complex formation of LD with water have been performed. 
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1. Introduction 



Levodopa (LD) is used as a pro-drug to increase dopamine levels for the treatment of Parkinson's 
disease [1]. Dopamine, itself cannot cross the blood-brain barrier but LD is able to. Once levodopa has 
entered the central nervous system, it is metabolized to dopamine by aromatic-L-amino-acid 
decarboxylase. LD is converted to dopamine in the body before reaching the brain. LD is metabolized 
by four major pathways: decarboxylation, O-methylation, transamination and oxidation as displayed in 
Chart 1 [2,3]. The principal path is decarboxylation, whereby dopamine is formed by aromatic amino 
acid decarboxylase [4-10], (cf. Chart 1). 

Chart 1. Two pathways of Levodopa (LD) metabolism. 



CH 3- 0 



OH 



OH 




OH 



(1) Decarboxylation 



OH 
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NH 2 



Numerous experimental methods have been used to analyze and determine the products of LD 
metabolism such as radioenzyme [11,12], and chromatographic [13] liquid chromatography (HPLC), 
UV [14], mass [15,16], electrochemical [17-20] and fluorescence [21-28]. The literature does not seem 
to contain any theoretical investigation of LD metabolism and as a consequence there is no systematic 
study of the metabolic reaction pathways, its energetics and transition state. 
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We have launched a research project with the ultimate aim of contributing to the understanding 
of the metabolic pathways of LD. This project involves a series of experimental and theoretical 
investigations of the different aspects of the molecular electronic structure of LD, acid-based 
properties and metabolic pathways of the decarboxylation processes of LD. 

In the present article, equilibrium geometry and charge density distribution of LD will be presented 
and discussed. Charge migration into or away from the side chain will be examined using natural orbital 
analysis. The important chemical and biochemical processes of protonation and deprotonation of LD 
will be examined. All possible protonation/deprotonation sites will be investigated; energetics and 
consequent charge redistribution will be computed, analyzed and discussed. The DFT method at a high 
level of theory will be used throughout this work. For a more quantitative understanding of the forces 
that govern the structure of LD, we will analyze the NBO orbitals and the second order perturbation 
energy stabilization due to internal donor-acceptor orbital interactions (principal derealizations). 

2. Results and Discussion 

2.1. Equilibrium Geometry of LD and Its Deprotonated Forms 

Figure 1 presents the geometry and numbering system of LD adopted in the present work. The 
DFT/B3LYP/6-311+G** equilibrium geometric parameters of LD and its deprotonated forms are 
presented in Figure 2. 

Figure 1. Numbering system of LD adopted in the present work and optimized geometric 
parameters computed at different levels of theory. 




Syn-MP2 (E opt = -703.567554 au) 
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1.462 




Syn-WB97XD (E opt = -705.453672 au) 




Syn-B3LYP (E opt = -705.45843 au) 
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Figure 1. Cont. 




Anti-B3LYP (£ opt = -705.453672 au) 



Figure 2. Geometrical parameters (bond lengths (A) and natural bond orbital (NBO) 
charge distribution) for deprotonated structures of levodopa. 



DP-1 
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DP-2 



DP-2-HB 



DP-3 
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Figure 2. Cont. 



DP-4 



DP-4 -HB 




Results of the present work, indicate clearly that LD is not planar, with the aliphatic side chain acting 
as a free rotator around several single bonds. The potential energy surface is broad and flat. This makes 
the process of arriving at the global minimum difficult and time consuming. Two shallow minima were 
identified to correspond to equilibrium ground state structures (cf. Figure 1); with a barrier height of 
few kcal/mol separating them. These two structures correspond to the syn and anti conformers 
resulting from the rotation of the OH group of the carboxyl group. The syn (LD-HB) conformation 
seems to be more stable by 2.987 kcal/mol due to the formation of H-bond with the amino nitrogen 
atom as proton acceptor. This H-bond stabilizes the syn conformation. 

The deprotonated forms of LD considered in the present work are DP-1, DP-2, DP-3, and DP-4 
(cf Figure 2). The geometries of all deprotonated forms have been fully optimized at the DFT/B3LYP/ 
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6-31 1+G** level of theory. Stationary points were characterized by frequency calculations at the same 
level of theory. Due to the low barriers to rotations around single bonds in the aliphatic side chain, the 
geometry optimization processes may very well correspond to one of the local minima other than the 
global minima. However, we have conducted molecular mechanics using the MM+ force field to 
explore the conformational space around the minimum. This search indicated that the potential energy 
surface is very shallow and the barriers are extremely low. Although a thorough exploration and analysis 
of the conformational space is required, we can still assume that the structures depicted in Figure 2 are 
very close to the global minima. It is interesting to note the effect of deprotonation on the general features 
of the equilibrium structure of LD. By inspection of Figure 2 the following remarks may be made: 

(1) In all cases studied, deprotonation is local in nature, i.e., its effect is localized in the region 
where deprotonation took place and is not transmitted through the length of the molecule. 

(2) Deprotonation of the carboxyl group (hydrogen atom H23) causes a derealization of charge in 
the 0-C=0 region. The two Rc-o bonds are considerably affected. Thus, the C-0 bond length 
is reduced from a value of 1.359 A to 1.256 A which is very close to the other RC-0 of 
1.251 A. The Cio-Cn bond length is slightly stretched by only 3% of its original value. This is 
accompanied by widening of the O12-C11-O13 angle to a value of 129.5° to minimize the 
repulsion due to the accumulation of the charge density in this region. The hybridization 
scheme of the oxygen atom remains, however, almost unaffected as: core 2s 1 ' 7 2p 5 ' 02 3p 001 

(3) It is of special interest to note that the effect of deprotonation of the carboxyl group proton is 
transmitted to the NH 2 group. The charge density on N is markedly increased. This point will 
be discussed further since it might very well underlay the formation of zwitterions. 

(4) In LD, the aliphatic side chain is tilted out of plane by an angle of 79°. Upon deprotonation of 
the carboxyl proton, the dihedral angle is reduced to 54° and the side chain is forced back to 
approach the plane of the rest of the molecule. 

(5) Deprotonation of the amino group has an even more localized effect. Thus, while the C-N bond 
length is slightly reduced, all other bond lengths are hardly affected. The charge is now 
accumulated on the N atom. This is not even transmitted to the C-N carbon. This is accompanied 
by a change in the N atom hybridization scheme from 2s L402 p 4413 P°' 01 to 2s L552 p 4 ' 393 P° 01 - The 
major effect is the accumulation of extra 1 1% e into the 2 s space. 

(6) Deprotonation of the catechol hydroxyl group causes the expected decrease in the phenyl-0 
bond length and increase of the negative charge on the oxygen atom. However, it seems that a 
subsequent hydrogen bond is formed with O-catechol hydroxyl H atom. This is stabilized by 
the increased negative charge on the O atom and the formation of a five-membered ring. The 
length of this H-bond is typical (1.881 A). 

2.2. Energetics of the Protonation/Deprotonation of ID 

There are three main sites of deprotonation (possibly 5 protons). The catechol hydroxyl group 
protons, the amino group protons and the carboxyl group proton. On the other hand, LD can abstract 
proton and act as base. The most possible site of protonation is the amino group nitrogen atom. In the 
present section we will investigate the proton affinity/proton detachment energies of LD. All possible 
sites will be examined. 
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Proton affinity, PA, or the deprotonation enthalpy, DE, of an acid in the gas phase, may be defined as 

A + H + = HA; PA = -AH° 29 8 

A small negative AH°298 value indicates strong acid. 
The enthalpy of deprotonation can be explicitly written as 

DE = -A E 0 elec + AZPE - | 0 298 ACp At 

Where A E 0 eXec is the difference in electronic energies of the acid and its deprotonated species, AZPE 
is the corresponding difference in zero point energies and J 0 ACp dt is corresponding difference of 
constant pressure heat capacities. Assuming that A~ and HA both possess the same number of degrees 
of freedom then the last term is reduced to the H + contribution of 5/2 kT (6.2 Kj/mol). This small 
contribution is omitted in tabulating the results of the present work. 

The relative energies of all deprotonated forms are included in Table 1 . Results of the present work, 
suggest that deprotonation of the NH 2 group is less favored than both the carboxylic group and catechol 
OH groups by 44.51 and 42.81 kcal/mol, respectively. It is interesting to note that the acidity behavior of 
these two latter hydroxyl groups is almost the same. As we can see in Figure 2, LD deprotonation to 
form DP-1 is predicted to be about 1.67 kcal/mol more stable than DP-4. These results indicated that 
the NH2 group is of the highest basicity. 



Table 1. Energetics of the Protonation/deprotonation of LD computed at the B3LYP/ 
6-3 11+G** level. 



Compounds 


EJnu 


Deprotonation and protonation 
energies (DE)/kcal/mol =£W - E e 


Relative energy 
kcal/mol 


LD 


-705.45367 




-2.984 


LD-HB * 




-705.45843 


0 


DP-1 


-704.90918 


-344.653 




DP-2 


-704.83826 


-389.156 


-43.496 




DP-2-HB * 


-704.90757 


-345.660 


DP-3 


-704.14807 


-822.248 




DP-4 


-704.88462 


-360.067 


-13.719 




DP-4-HB * 


-704.90648 


-346.348 


P07 


-705.74793 


181.6252 




PN 


-705.82031 


227.0671 




P013 


-705.81224 


221.9833 





* HB means existence of intramolecular hydrogen bond. 



It should be noted that protonation/deprotonation energies reported here are corrected for electron 
correlation. Even though this point is very important, it is ignored in most PA/PDE computations. By its 
very definition, reactions under investigations involve proton transfer where the total number of electron 
spins does not change, i.e., the reaction is isogyric proton transfer. For such a reaction, correlation effects 
are very important. 

It is interesting to carefully examine the deprotonation of carboxylic group and that of the amino 
group, DP-1 and DP-2. In case of DP-1 deprotonation took place in one step followed by derealization 
of charge in the 0-C=0 region. NH 2 group deprotonation, however, involves three distinct steps. First, 
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the deprotonation of the amino group proton followed by rotation of the OH group around the C-0 bond 
and finally, transfer of the hydroxyl group proton to the nitrogen atom. H-bond formation stabilizes this 
final structure considerably (E = 49.5 kcal/mol). Deprotonation of the aromatic ring hydroxyl group 
leads to a structure that is stabilized by moderate intensity H-bond (13 kcal/mol). 

In conclusion, deprotonation reactions considered in the present work seems to involve two 
subsequent steps. First, deprotonation, which is followed by rearrangement to form H-bonded structures. 
This is evident in the case of DP-2 and DP-4. Deprotonated forms seem to enjoy extra stability due to 
intramolecular H-bonding (Figure 2). Three possible sites of protonation of LD were considered in the 
present work, the ring hydroxyl oxygen atom, the side chain amino nitrogen and oxygen atom of 
carbonyl group. Figure 3 presents the optimized structure of the three protonated forms. Results of the 
present work, indicates clearly that protonation of the amino acid side chain is more favorable. It is 
interesting to consider in detail the protonation of the amino acid hydroxyl oxygen. As Figure 3 shows, 
protonation of the oxygen atom is followed by transfer of a proton to the nitrogen of the amino group 
this is followed by 1,3-proton shift in the carboxyl group and derealization of the 7i-charge density in 
the 0-C=0 region. This suggests that the nitrogen atom in LD acts as strong proton acceptor, a very 
crucial step in zwitterions formation. 

Figure 3. Geometrical parameters (bond lengths (A)) of the protonated forms for LD 
computed at the B3LYP/6-31 1+G** level. 



P07 
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Figure 3. Cont. 




2.3. Water-LD Complexes 

The B3LYP/6-311+G** optimized geometries of the three cyclic LD-water complexes DL-W1, 
DL-W2 and DL-W3 are shown in Figure 4. 

Let us first check the effect of hydrogen bonding on the bond length in monomers. In the case of 
LD-W1, it can be seen that the distance C11-O13 is increased, whereas C11-O12 is decreased 
compared to the corresponding distance in Levodopa monomer. On the other hand, the distance 
C11-O13 is decreased, whereas C11-O12 is increased in LD-W2 form. There are slight changes in the 
other bond lengths in all three water complexes. 
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Figure 4. Optimized structures for the LD- water (distances in A). 




Table 2 presents the total energies (£e), zero point energies (ZPE), binding energies A£b and corrected 
binding energies Afi^ZPE) of LD-water complexes computed at the DFT/B3LYP/6-311+G** level. 
Table 2 also includes the energy corrected by the counter poise method for the BSSE. Inclusion of this 
energy term changes the order of stability of the LD-W complexes. It is interesting to notice that not 
all LD-W complexes suffer from the BSSE to the same extent. Thus, LD-W1 shows a maximum BSSE 
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of 10.44 kcal/mol, whereas LD-W2 shows a minimum BSSE of 3.96 kcal/mol. Results reported in 
Table 2 reveal that the LD-W2 complex is slightly more stable than LD-W1. The most stable complex 
(LD-W3) is formed at the site characterized by the highest acidity (OH). The strongest bonds occur 
when hydrogen bridges form between oxygen atoms acting as donor and acceptor. The hydrogen bond 
undergoes elongation by ca. 0.36 A when the amino group becomes the proton donor. We can explain 
the stability of these complexes in terms of the size of the ring formed. In the case of the 
OH group acting as a proton donor, the hydrogen bond is a part of a six-membered ring (Figure 4), its 
enhanced strength is due to reduction of the ring strain. In contrast, when the amine group participates, 
as H-donor, and the OH group as acceptor, in the hydrogen bond formation, a seven-membered ring is 
formed (Figure 4). This ring contains long and week hydrogen bonds. 



Table 2. Total energies, enthalpies and relative energies of levodopa-water complexes 
computed at the B3LYP/6-31 1+G** level. 



Compounds 


£ e /au 


(AE h ) 
kcal/mol 


^BSSE 


(A^b-bsse) 

kcal/mol 




(A#eb) 

kcal/mol 


7/\au 


AH 
kcal/mol 


L-Dopa 


-705.45367 








-705.25677 




-705.24186 




LD-W1 


-781.92896 


-10.559 


-781.91194 


-0.11922 


-781.70706 


-8.206 


-781.68921 


-8.788 


LD-W2 


-781.91946 


-4.598 


-781.91112 


-0.63377 


-781.69836 


-2.783 


-781.68007 


-3.026 


LD-W3 


-781.92633 


-8.905 


-781.90779 


-2.72335 


-781.70373 


-5.939 


-781.68652 


-7.073 


Water 


-76.45846 








-76.43716 




-76.43339 





E Q : Energy of optimized structure; AE h : Binding energy; E c = E Q + ZPE AE ch : Corrected binding energy; 



H: Enthalpy; ^bsse- counter poise corrected energy. 

The binding energies for LD-water complexes follow the order: LD-W1 > LD-W3 > LD-W2, 
which indicates that when the hydroxyl groups act as H-donor, the binding energy of the complex will 
increase as compared to the case when the NH2 is the H-donor. Data in Table 2 suggest that the energy 
barrier for the decarboxylation process may very well decrease upon interaction of LD with water. 

2.4. NBO Analysis 

The natural bond orbital analysis provides an efficient method for studying intra- and 
inter-molecular bonding and interaction among bonds, and also provides a convenient basis for 
investigating charge transfer or conjugative interaction in molecular systems. The interactions between 
"filled" (donor) Lewis-type NBOs and "empty" (acceptor) non-Lewis NBOs lead to loss of occupancy 
from the localized NBOs of the idealized Lewis structure into the empty non-Lewis orbitals, and they 
are referred to as "derealization" corrections to the zeroth-order natural Lewis structure. The 
stabilization energy AEy (kcal/mol) associated with derealization is estimated by the second-order 
perturbative as 

AEg = qi(F (iJ) 2 )/(sj - St) 

where q t is the donor orbital occupancy, e,-, ey are diagonal elements (orbital energies) and Fqj) is the 
off-diagonal NBO Fock matrix element. 



Int. J. Mol Sci. 2012, 13 



4334 



NBO analysis has been performed on LD at the DFT/B3LYP/6-31 1+G** level in order to elucidate 
the derealization of electron density within the molecule, and is presented in Table 3. The stabilizing 
and destabilizing energies of NBO interactions were considered up to 1.2 kcal/mol. 



Table 3. The second-order perturbation energy (derealization) calculated for LD and 
its anions at B3LYP/6-31 1+G** level. 



Compounds 


Donor 


Acceptor 


E (2) 




LP 07 


7I*C1-C2 


26.47 




LP 08 


ft*C3-C6 


23.54 


L-Dopa 


LP 012 


TC*C11-013 


41.99 




LP 013 


rc*ci 1-012 


34.54 




LPN14 


G*C10-C11 


8.80 




LP 07 


rc*ci-c2 


23.95 




LP 08 


7I*C3-C6 


20.03 


DP-1 


LP 012 
LP 013 


rc*ci 1-013 

TC*C11-012 
G*C10-C11 


18.55 
19.45 
18.32 




LPN14 


TC*C9-C10 


3.63 




LP 07 


7t*Cl-C2 


23.07 




LP 08 


7I*C3-C6 


19.27 




LP 012 


JI*C1 1-013 


28.99 


DP-2 


lp on 


7I*C11-012 
<7*C10-C11 


33.59 
15.74 




LPN14 


TC*C11-013 
<7*C10-C11 


8.31 
10.78 




LP 07 


7I*C1-C2 


21.90 




LP 08 


ft*C3-C6 


12.90 




LP 012 


t*Cl 1-013 


17.88 


DP-3 


LP 013 


ft*Cll-012 
t*C10-Cll 


17.42 
18.50 




LPN14 


O*C10-H22 
G*C10-C11 


10.04 
8.63 




LP 07 


ft*Cl-C2 


27.89 




LP 08 


7l* C 3-C6 


71.04 


DP-4 


LP 012 
LP 013 


ft*Cll-013 
ft*Cl 1-012 
O*C10-Cll 


38.99 
34.28 
16.05 




LPN14 


G*C10-C11 


8.96 



Carful inspection of the NBO energy analysis of LD reveals the following: 

1 . The major D-A interaction in LD is localized on the phenyl ring and involves 7rc-c -7l *c-c 

2. Strong interaction between the p lone-pair electrons of the hydroxy 1 group oxygen atoms and 
the ring 7r*c-c- This CT interaction spans a narrow range of 23-26 kcal/mol. 
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3. A very strong interaction involving the p lone pair of the carboxyl group O12 and the 7C*ci 1-013. 
This interaction amounts to 41.99 kcal/mol and reflects the strong derealization of the n- electron 
density in the O-C-0 region. 

4. The p lone pair of O13 is, surprisingly interacting with the Rydberg orbitals of the carbonyl 
carbon. This lone pair also interacts via CT to 7i*c-oi2 by a considerable amount of energy of 
34.54 kcal/mol. 

5. In conclusion, LD behaves as if it is composed of two non- or at least weekly interacting 
subsystems, namely the aryl moiety and the amino acid side chain. One can trace, however, 
week interaction through the a- framework of the aryl moiety transmitted to Oc9-cio which 
interacts with g*cio-ni4 and to a lesser extent with 0*011-013. Furthermore, Oc9-cio also interacts 
with the Rydberg type orbitals of Cn and N14. Thus, although, there is no n conjugation 
extending over the LD framework, yet it seems that the a framework is able to transmit the 
interaction between the catechol moiety and the amino acid side chain. 

Now, it is interesting to examine the effect of deprotonation on the charge density distribution and 
bonding interactions in LD. Inspection of NBO analysis of the deprotonated forms reveals the following: 

1. DP-1, where the carboxyl group proton has been detached, shows a dimensioned interaction 
across the molecular framework. Thus, the interaction of Oc9-cio with the amino group side 
chain is reduced to almost half its value in LD. However, it maintains its week interaction with 
the Rydberge type orbitals on Cn and N14. Furthermore, a new very strong interaction appears 
between each of the lone pair orbitals on O12 and O13 with the empty antibonding p orbital on Cn. 
This interaction amounts to 191 and 205 kcal/mol for O12 and O13 respectively. This indicates 
clearly that deprotonation of the carboxyl group has a pronounced effect on the charge density 
distribution that is localized to a great extent in the 0-C=0 region. 

2. For DP-2 where the amino group proton is detached the process has very little effect on the 
charge density distribution. Thus, the across- subsystem interactions through Gc9-cio is enhanced 
slightly. A new week interaction emerges between the N14 lone pair and the antibonding g*cio-cii 
orbital. Deprotonation in this case has the subsequent effect of localizing the charge density onto 
the nitrogen atom which acts as a very poor electron donor. 

3. For DP-4, where the catechol moeity hydroxyl proton is removed the deprotonation process has 
strong perturbation effect on the 7i charge density of the aryl moiety. Thus, a strong interaction 
appears between the deprotonated O13 lone pair and the 7i*c3-c6 This interaction (E = 71 kcal/mol) 
is delocalized over the entire catechol n- framework, and has no marked effect on the 
interaction between the aryl moiety and the amino acid side chain. 

NBO analysis provides a quantitative predictive tool for the protonation/deprotonation of LD. Thus, 
let us analyze the energetics and the population of the lone pair orbitals on oxygen and nitrogen 
(cf. Table 4). The four oxygen atoms, in LD, each have two lone pair orbitals, one of lower energy of 
the sp hybrid type and the other is a pure p orbital. The carbonyl group oxygen has the highest energy 
lone pair p type orbital which thus represents the most probable site of protonation. The amino group 
nitrogen carries one lone pair which is a hybrid sp type orbital. It should be noted that this LP is of 
over 80% p character, the s-contribution would lead to greater symmetry and extension in space. This 
most probably underlies our previous observation that nitrogen is a better proton acceptor. 
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Table 4. NBO energy and hybridization schemes of the O and N lone pairs and O-H and 
N-H bond orbitals in LD. 



Occupancy 


Type 


Atom 


Hybridization coefficients 


Energy, au 


(1.97866) 


LP(1) 


07 


s(45.05%)p 1.22(54.91%) 


-0.60588 


(1.87642) 


LP (2) 


07 


s(0.00%)p 1.00(99.94%) 


-0.32010 


(1.97737) 


LP(1) 


012 


s(44.86%)p 1.23(55.11%) 


-0.64468 


(1.82991) 


LP (2) 


012 


s(0.00%)p 1.00(99.94%) 


-0.35517 


(1.97844) 


LP(1) 


013 


s(58.78%)p 0.70(41.20%) 


-0.71385 


(1.94184) 


LP(1) 


N14 


s(18.42%)p 4.43(81.53%) 


-0.31638 


(1.98704) 


BD(1) 


07-H18 


0(74.69%) H(25.31%) 


-0.73784 


(1.98734) 


BD(1) 


012-H23 


0(74.86%) H(25.14%) 


-0.77110 


(1.98866) 


BD(1) 


N14-H24 


N(68.12%) H(31.88%) 


-0.62658 



Table 4 presents the NBO results of the O-H and N-H bond orbitals. O-H bonds seem to be more 
ionic in character which would facilitate deprotonation of the O-H bonds. It should be noted that, 
Table 4) presents non = redundant enters only, that is one O lone pair, one O-H bond pair and so forth. 

In conclusion, NBO analysis provides a powerful predictive tool for the protonation/deprotonation 
processes. The analysis of this tool depends on the energetic, population and extent of hybridization of 
the LP or X-H bond orbitals. 

3. Computational Methods 

All calculations have been be carried out using the Gaussian03 [29] package of programs. The 
geometries of LD and its anions have been fully optimized at the DFT/B3LYP/6-311+G** level of 
theory [30-32]. Frequency calculations were performed at the same level of the theory in order 
to characterize the stationary points and to evaluate the zero-point energy (ZPE). The stable LD 
conformation was subjected to further confirmation and validation using another functional, namely 
WB97XD and MP2 method with the same basis set. Results of the present work indicate clearly that 
the B3LYP is more reliable than the WB97XD and MP2. Thus, all calculations carried out in the 
present work are at the B3LYP/6-31 1+G** level of theory. The LD-water complexes were constructed 
starting from the most stable LD conformer. All complexes were also fully optimized at the same level, 
B3LYP/6-311+G**. Interaction energies in LD-water complexes may be affected by the basis set 
superposition error (BSSE). The counterpoise (CP) method of Boys and Bernardi [33] has been used to 
correct for the BSSE. The natural bonding orbitals (NBO) calculations [34] were performed using 
NBO 3.1 program as implemented in the Gaussian 03 package at the DFT/B3LYP/6-31 1+G** level in 
order to understand various second-order interactions between the filled orbitals of one subsystem and 
the vacant orbitals of another subsystem, which is a measure of the inter-molecular derealization or 
conjugation. 

4. Conclusion 

The present work represents the foundation for a thorough and detailed investigation of the 
mechanism of decarboxylation of levodopa, a drug of vital importance which is used to increase 
dopmaine level for treating Parkinson's disease. The major metabolic step of LD to produce dopmaine 
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is the decarboxylation. In the present work we report the equilibrium geometries of LD and all its 
possible protonated and deprotonated forms computed at the DFT/B3LYB/6-31 1+G** level of theory. 
Results of the present work indicate that: 

1. LD is not planar with its side chain acting as a free rotor across several single bonds. However, 
deprotonation of the carboxyl group or the amino group forces the side chain into the plane of 
the catechol moiety. 

2. Deprotonation is local in nature, for all sites considered the changes in geometry and charge 
density are localized in the region where the deprotonation took place. However, the effect of 
carboxyl group deprotonation is transmitted to the C-NH 2 region and leads to the accumulation 
of excess charge density on N14. This point is of special importance when considering the 
formation of zwitterions. 

3. Deprotonation of the amino group proton lead to the change of the hybridization scheme of N14. 

4. NBO analysis of LD and its deprotonated forms reveals detailed information of bonding 
characteristics and interactions across the molecular frame work. Thus, although there is a 
pronounced cross conjugation in the 7i-system, the interaction between the catechol moiety and 
the amino acid side chain is marked through the inductive and CT effects across the a and 
nonbonding frameworks. 

5. Deprotonation of the carboxyl group is more favorable than deprotonation at other sites. 
However, it should be noted that stabilization of the deprotonated structures by intramolecular 
H-bonding is visible in the case of the deprotonation of the amino group or the catechol 
hydroxyl group. 
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